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ABSTRACT 

Young, rapidly rotating neutron stars could accelerate ions from their surfaces to 
energies of ~ 1 PeV. If protons reach such energies, they will produce pions (with 
low probability) through resonant scattering with x-rays from the stellar surface. The 
pions subsequently decay to produce muon neutrinos. Here we calculate the energy 
spectrum of muon neutrinos, and estimate the event rates at Earth. The spectrum 
consists of a sharp rise at ~ 50 TeV, corresponding to the onset of the resonance, 
above which the flux drops with neutrino energy as e~^ up to an upper-energy cut-off 
that is determined by either kinematics or by the maximum energy to which protons 
are accelerated. We estimate event rates as high as 10-100 km~^ yr""'^ from some 
candidates, a flux that would be easily detected by IceCube. Lack of detection would 
allow constraints on the energetics of the poorly-understood pulsar magnetosphere. 
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1 INTRODUCTION 

Astrophysical neutrinos of high energy ( <; 1 GeV) are ex- 
pected to arise in many environments in which protons 
are accelerated to relativistic energies. Neutrinos produced 
by the decay of pions created through hadronic interac- 
tions (pp) or photomeson production (^17) escape from 
the source and travel unimpeded to Earth, and carry in- 
formation directly from the acceleration site. Neutrinos 
may be produced by cosmic accelerators, like those i n su- 
pernova remnants iProthBroej_^edrwek_fc Luo[ Jl998l) . ac- 
tive gala ctic nuclei J Learned fc MannheimI l2000f) . micro- 
quasars JDistefanO-et_aL 2002i) and ga mma-rav bursts 
JWaxman fc BahcalllllOOTl : iDai fc LulbOOlTl . To detect these 
neutrinos, several projects are underway to develop large- 
scale neutrino detectors under water or ice; AMA NDA-II, 
ANT ARES and Baikal are running, while IccCube JHalzenl 
l2006t) . NEMO and NESTOR (lCarin2003) are under con- 
struction. 

Recently, we proposed that young {tags < lO^yr) and 
rapidly-rotating neutron stars could be intense neutrino 
sources (Link fc Burgio 2005; hereafter LB05). For stars 
that have a stellar magnetic moment with a component anti- 
parallel to the spin axis (as we expect in half of neutron 
stars), ions will be accelerated off of the surface; otherwise, 
electrons will be accelerated. If energies of ~ 1 PeV per pro- 
ton are attained, pions will be produced through photome- 
son production as the protons scatter with surface x-rays. 
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producing a beam of ^ neutrinos with energies above ~ 50 
TeV. Detection of such neutrinos would provide an invalu- 
able probe of the particle acceleration processes that take 
place in the lower magnetosphere of a neutron star. In this 
paper we predict the neutrino spectrum that would result 
from this production mechanism to aid in the interpreta- 
tion of experimental results from searches for astrophysical 
neutrinos. We obtain improved estimates of the event rate. 
In the next section, we review the acceleration model 
of LB05. In Section III, we calculate the neutrino spectrum 
that results from this model. In Section IV, we estimate 
the count rates that could be seen in a km^-scale experi- 
ment such as IceCube. We conclude with a discussion of the 
prospects for detection of high-energy neutrinos from pul- 
sars. 



2 THE MODEL 

In the neutrino production scenario of LB05, protons (within 
or without nuclei) are accelerated in the neutron star magne- 
tosphere to high enough energies to undergo resonant scat- 
tering with surface x-ray photons (the A resonance): 
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The A+ is a short-lived excited state of the proton with a 
mass of 1232 MeV. For this process to be effective, ions must 
be accelerated close to the stellar surface, where the photon 
density is high and the process is kinematically allowed (see 
below, and LB05) . The plasma is tied to the magnetic field, 
so the acceleration can occur only in the direction of the 
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magnetic field B. In a quasi-static magnetosphere with a 
magnetic axis that is parallel or anti-parallel to the rotation 
axis, the potential drop across the field lines of a star ro- 
tating at angular velocity Q — 2-k jp (where p is the period) 
from the magnetic pole to the last field line that opens to 
infinity is of magnitude JGoldreich fc Julianlll969r : 



A*^ 



2c2 



7 X 10^**Bl2i?6Pms Volts, 



(2) 



per ion of mass number A and charge Z. Here B = 10^^i3i2G 
is the strength of the dipole component of the field at 
the magnetic poles, R = IO^-Rb cm is the stellar ra- 
dius and pnis is the spin period in milliseconds. Hence- 
forth we take R(, = \ when making estimates. In equi- 
librium (not realized in a pulsar), a co-rotating magneto- 
sphere would exist in the regions above the star in which 
magnetic field lines close; the charge density would be 
pq ~ e^no — Bjpc (cgs), where no is the Goldreich- 
Julian number density of ions. Deviation from corota- 
tion will lead to charge-depleted gaps somewhere above 
the stellar surface, through which charges will be acceler- 
ated to relativistic energies llRuderman fc Siitherlandlll975t 
lArons fc Scharlc mannlll979r . The proton energy threshold 
tp for A""" production is given by 



cos^)"^ 0.3 GeV^ 



(3) 



where e^, is the photon energy and Q is the incidence angle 
between the proton and the photon in the lab frame. Young 
neutron stars typically have temperatures of Too — 0.1 keV, 
and photon energies e^ = 2.%kTrx(\ -f Zg) ~ 0.4 keV, where 
Zg ~ 0.4 is the gravitational red-shift (for i? = 9 km and 
M — I.4M0) and Too is the surface temperature measured 
at infinity. For 6 ~ n/2, corresponding to scattering of pro- 
tons near the stellar surface with photons from the stellar 
horizon, the proton threshold energy for the A resonance is 
then ep,th ==; r-j^kcv PeV, where To.ikoV = (fcToo/O.l keV). 
Let us compare the required energy for resonance to the 
potential drop per proton across B: 
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Hence, for a pulsar spinning at 10 ms, a potential drop along 
B of only ~ 1% of that across B will be sufficient to bring 
protons or low-mass nuclei to the resonance. Optimistically 
assuming that the full potential A"!? is available for accel- 
eration along field lines, a necessary condition for the A 
resonance to be reached is {A = 1) 

-Bi2P,;Jro.ikeV ^ 3 X 10-^ (5) 

Assuming To.ikeV ~ 1, typical of pulsars younger than ~ 10^ 
yr, there are 10 known pulsars within a distance of 8 kpc 
that satisfy this condition (Manchester et al. 2005), about 
half of which should have positively-charged magnetic poles; 
these are potentially detectable sources of ^ neutrinos. The 
best candidates are young neutron stars, which are usually 
rapidly spinning and hot. Equality corresponds to the full 
potential across field lines being present along field lines, 
probably an unlikely scenario, since space charge will act 
to quench the electric field along B. For stars in which the 
above inequality is easily satisfied, the protons will reach 
energies sufficient to undergo photomeson production if the 
electric field along a typical open field line is much smaller 
than the field across the line, the more probable situation. 



In the photomeson production process of eq. Q, each 
muon neutrino receives 5% of the energy of the proton. [We 
assume, for simplicity, that the pious do not undergo sub- 
sequent acceleration]. Typical proton energies required to 
reach resonance are ~ TfTikcV P^V, so the expected fj, neu- 
trino energies will be e,, ~ ^^'^omev TeV. Moreover, since 
the accelerated protons are far more energetic than the ra- 
diation field with which they interact, any pious produced 
through the A resonance, and hence, any muon neutrinos, 
will be moving in nearly the same direction of the protons. 
The radio and neutrino beams should be roughly coincident, 
so that some radio pulsars might also be detected as neu- 
trino sources. We see the radio beam for only a fraction ft 
of the pulse period, the duty cycle. Typically, /;, ~ 0.1 — 0.3 
for younger pulsars. We take the duty cycle of the neutrino 
beam to be ft (but see below). In LB05, we estimated the 
phase-averaged neutrino flux at Earth resulting from the ac- 
celeration of positive ions, at a distance d from the source, 
to be 
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where fd is the fraction by which the space charge in the 
acceleration region is depleted below the corotation density 
no, and Pc is the conversion probability. A better estimate 

is 
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(7) 



The pre- factor /ti(l — /d) is a more realistic interpolation 
between the two regimes of complete depletion {fd ~ 0) 
and no depletion {fd = 1). In the latter case, there should 
be no neutrino production as the field along B is entirely 
quenched. The factor of two arises from inclusion of the t-p 
from eq. Q. For the purposes of calculating the spectrum, 
we use the following differential form for the neutrino flux: 



— ^2cf,fd{l^fd)noi- — . 



(8) 



Here we are interested in estimating upper limits on the flux; 
we henceforth take fd = 1/2 and Z — A — 1. 

In LB05, we estimated Pc by assuming that the protons 
reach the energy required for photomeson production imme- 
diately above the stellar surface; we then evaluated Pc at a 
particular height above the stellar surface. The focus of this 
paper is to obtain the spectrum dPc/dev We account for 
the fact that the proton acceleration will take place over a 
finite distance above the surface. We also include the finite 
width of the A resonance in the cross section. 



3 NEUTRINO SPECTRUM 

Let the accelerated proton have energy e^ and the sur- 
face photons have energy e-y. Define a dimensionless energy 
X — epe.y/eo, where eq = 0.3 GeV^. To account for the fi- 
nite width of the A resonance, let us express the resonant 
contribution to the energy-dependent p7 cross section for a 
proton of energy e in the nucleon rest frame as 



o" = (To cxp 



2«)2 



for e > eT ^ w, 



(9) 
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Figure 1. Geometrical setup. See text for details. 



where cro = 5 x 10~ cm~ is the cross section for A"*" 
production, ey is the A energy threshold and w ~ 100 MeV 
is the width of the A resonance. For a given photon energy 
and scattering angle, the proton energy at the resonance is 
given by 



ej, = IO(l„cos0)-' 



XT 



(1 



(10) 



The cross-section is non-zero for x > xt — x^ where x^ = 
we^/eo — 0.1. This implies a minimum kinematic scattering 
angle 9o given by 



cos 6o = 1 — {x + Xw) 



(11) 



In Fig. 1 we show the geometry that we use for calculating 
the spectrum. We neglect the effects of gravitational light 
bending, which would act to increase the rates we calculate 
here. The proton is located at the point denoted by p, at 
a height z above the stellar surface, and 9 is the scattering 
angle between the proton and the incoming photon. The 
scattering angle 6 is always less than 7r/2, which requires 
the protons to have energies that satisfy a; > 1 — a;^, in order 
to undergo resonant scattering. There is also a maximum 
scattering angle 6h defined by the stellar horizon as seen 
from height z: 



sm9h = {1 + z) , 



(12) 



where z is the height above the stellar surface in units of 
the stellar radius R. We must have 9o < 9^ for resonance to 
occur. We can rewrite the cross-section as 



o"(a;, 9) = o"o exp 



-(a; - xtY 



for X > xt — Xu 



(13) 



To calculate the conversion probability, we need the number 
density of photons as a function of height z and incident 
photon angle 6. Let the number density of photons at the 
stellar surface be n-y(O). An element of surface radiates into 
2-n steradians, independent of angle. For a small surface el- 



ement of area dA, the contribution dn to the total photon 
number density n is 



n..{Q)dA . , 
dn = TT n ■ r 



h ■ f = cos(a + 6). 



(14) 



27r r2 

Taking the ring surface element shown in the figure, that 
area is 

dA^2-KR^smada. (15) 

Some useful trigonometric identities are 

rsin& = _Rsina (l+z) sin6' = sin/? 13 = -n— a— 6, (Id) 

which give 

a = sm-^[{l + z)sm9]-9 (17) 

da (1 + z) cos9 



d9 v'l-(l + ^)'sm'6l 
The derivative has a singularity at 
-1 1 



-Ih = sm 



1 + Z-- 



(18) 



(19) 



which corresponds to the horizon angle of eq. 1121 . The num- 
ber density of photons at height z, arriving from angles in 
the range to 9 + d9 is 



, / „^ , ^ ,„ da sin^ 9 
dn{z,9) = n-y{0)d9-—— cos(a + 6*) 



(20) 



d9 sin a 

Suppose the protons have an energy that depends on the 
height over which they have been accelerated: e{z). At height 
z, the mean-free-path l{x, 9) for p —* A"*" conversion through 
scattering with a photon arriving at an angle between 9 and 
9 ^ d9 is given by 

r^{z,9)=dn{z,9)G{e{z),9). (21) 

To obtain the total mean-free-path for conversion at height 
z, we integrate over all possible angles for the incident pho- 
tons: 



r\z)=l\9'^.ieiz),9). 



(22) 



The probability for p — * A conversion of a proton as it 
moves from 2: to 2 -I- dz is 



dPc = TT^dz. 
l{z) 



(23) 



Let us now assume a specific dependence of the proton en- 
ergy on height: 

'')\ (24) 

where L is the characteristic acceleration length of the pro- 
ton in units of R. Then 



dz = —X 

7 



/7-1 



dx. 



(25) 



The probability of conversion per unit (dimensionless) en- 
ergy interval is 



dP, RLx^'-'-^ ["^ dn{x,9) . ^, 



(26) 



This is the energy spectrum of converted A^ particles. The 
neutrino spectrum is the same, but is scaled down in energy 
according to Xv = 0.05a::. The amplitude increases by a factor 
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Figure 2. The neutrino energy flux is displayed for the Crab and the Vela pulsars, for the cases of linear (solid line) and quadratic 
proton acceleration (dashed line). 



of 20 to maintain total probability. From eq. we find (in 



units of GeV m 



dev 



1.1 X 10-VoBi2P;;sdk' "^^^ 



kpc 



det. 



(27) 



where no = 7 x 10^^Bi2Pms cm ■^. We obtain dPc/dey from 
eq. 1261 . Since x,j = 0.05a;, we have 



dPc _ t~i dPc 



dtu 



dx 



dP 

2.7 X lO-^Toikcv-r^ GeV~ 
dx 



(28) 



Combining eq. 
flux: 



with eq. Q gives the neutrino energy 



— = 3 X 10 ^fbBl2Pmldy,p^To 



dP, 



de, ..^..^ms-Kpc- — V ^^ (29) 

Eq. 12911 is our chief result. It allows comparison with ob- 
served event distributions as well as an estimate of the total 
expected counts. In Fig.|21we show the neutrino energy flux 
for two candidate sources, the Crab and Vela pulsars, located 
respectively in the northern and southern hemispheres. For 
illustration, we consider linear (7 — 1) and a quadratic 
(7 = 2) proton acceleration laws. Linear acceleration cor- 
responds to an accelerating field that is constant in space. 
Quadratic acceleration corresponds to an accelerating field 
that grows linearly with height above the star. We have used 
x^o — 0.1. 

For either acceleration law, the spectrum begins sharply 
at 



70T,: 



TeV, 



(30) 



corresponding to the onset of the resonance. At higher ener- 
gies, the spectrum drops approximately as e^^, as the phase 
space for conversion becomes restricted; higher energy neu- 
trinos are produced by protons that have been accelerated 
to greater heights, where the photon density is lower and the 
solid angle subtended by the star (as seen by the proton) is 
smaller. At some maximum energy, the spectrum is suddenly 
truncated by either kinematics (solid curve) or the termina- 
tion of the proton acceleration as limited by the magnitude 
of the acceleration gap (not shown, since this cut-off has not 
been predicted). 

In obtaining eq. 1291 . we neglected the effects of gen- 
eral relativity, except when relating the stellar temperature 
at infinity to the temperature near the surface. General rela- 
tivistic effects bend the photon trajectories and bring some 



photons from beyond the classical stellar horizon (eq. 1121 
to where they can interact with the protons. This effect en- 
hances the rate, but not significantly. To estimate the mag- 
nitude of this effect, we can regard the star as effectively 
larger, and replace the stellar radius that appears as a pre- 
factor in eq. 1261 by the effective stellar radius at infinity, 
(1 -I- Zg)R. The effects of gravitational light bending will in- 
crease the flux by a factor of ~ 1 + Zg ~ 1.4 for R — 9 
km and a factor of only 1.2 for i? = 12 km. In light of the 
large uncertainties in the values of fd and L, we simply took 
R = 10 km in eq. 12(ill for the purpose of making estimates. 



4 ESTIMATED COUNT RATES 

We now use the spectrum obtained in the previous section to 
estimate the count rate in a detector. Large-area neutrino 
detectors use the Earth as a medium for conversion of a 
muon neutrino to a muon, which then produces Cerenkov 
light in the detector. The conversion probability in the Earth 
is (Gaisser et al. 1995): 



p 



1.3 X 10" 



_\ 



.1 TeV/ 
The muon event rate is 



dN 
dAdt 






(31) 



(32) 



The choice of the upper limit of integration is not very im- 
portant, because the spectrum is steep. Estimated count 
rates are given in the last column of Table 1 for the Crab, 
Vela and 9 other pulsars, assuming L = 0.1 and linear ac- 
celeration. Estimated conversion probabilities for p — > A"*", 
obtained by integrating eq. 12611 from xt to 10a;T, are shown 
in the penultimate column. Our event rates are a factor of 
~ 10 — 30 lower than estimated in LB05. The main reason 
for the lowered rate is phase space limitations imposed by 
the geometry; protons at a given height can only undergo 
resonant conversion from photons arriving from the surface 
in a narrow range of angles. Moreover, there is a competi- 
tion between the finite distance over which protons reach 
sufficient energy to be converted, and the reduction in the 
number of photons arriving from the star with the correct 
angular range for resonant conversion to be kinematically 
allowed. 



Source 



Crab 
Vela 

J0205+64 
B1509-58 
B1706-44 
B1823-13 
Cass A 
SN 1987a 



^kpc 



age 



Pins 



Bi2 To. 
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dN/dAdt 



2 

0.29 

3.2 

4.4 

1.8 

4.1 

3.5 

50 



10'^ 
1042 

102.9 

10^-2 

104.3 
104.3 

300 

17 



33 3.8 < 1.7 (Weisskopf et al., 2004) 0.14 

89 3.4 0.6 (Pavlov et al., 2001) 0.04 

65 3.8 ^ 0.9 (Slane et al., 2002) 0.05 

151 15 1? 0.26 

102 3.1 1? 0.13 

101 2.8 1? 0.34 

10? 1? 4 (Pavlov ct al., 2004) 0.1? 

1? 1? 4? 0.1? 



1.6 X lO-'^' 
7.2 X 10-5 
2.4 X lO-"* 
3.4 X 10-"' 
3.4 X lO-'' 
3.4 X 10-'' 
2.1 X 10-2 
2.1 X 10-2 



45 

25 

1 

5 

5 

2 

50 

3 



Table 1. Estimated upper limits on the /x fluxes at Earth. Numbers followed by question marks indicate guesses. Radio pulsar spin 
parameters were taken from the catalogue of the Parkes Radiopulsar Survey (Manchester et al. 2005). Temperatures and limits on 
temperatures were taken from the references indicated. The temperature upper limits on the Crab and J0205+64 were used. The 
integrated conversion probability Pc is reported for the case of linear acceleration, and assuming L = 0.1. 
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Figure 3. The muon event rate is shown as a function of the 
acceleration length for the Crab (circles) and the Vela (squares) 
pulsars. Upper (lower) curves display results for the case of linear 
(quadratic) proton acceleration. Solid and dashed lines guide the 
eye. 



In Fig. |3| we show the muon event rates estimated for 
the Crab (circles) and the Vela pulsars (squares), as a func- 
tion of the acceleration length L. The upper curves refer to 
calculations performed for a linear acceleration law, whereas 
the lower ones assume a quadratic one. The event rate in- 
creases up to a maximum value, and then decreases. Increas- 
ing L initially enhances the rate because the protons then 
attain enough energy to undergo resonant scattering at a 
height where much of the stellar surface is visible. If L is 
made too large, however, the process is reduced by a lower 
photon density and restricted range in photon angles avail- 
able for scattering. The dependence on L is more gradual for 
the quadratic acceleration model, because the protons must 
travel farther to attain enough energy for resonance. 



5 DISCUSSION 

To summarize, if protons reach the photomeson production 
resonance in the neutron star magnetosphere, they will pro- 
duce a spectrum of muon neutrinos with the following simple 
characteristics: 

(i) A sharp turn on at 



corresponding to the onset of the resonance. 

(ii) A rapid fall with energy as €~^ , determined by scat- 
tering kinematics. 

Neutrinos are produced at relatively high rates only 
if the protons are accelerated through the resonance close 
to the star {L'^ 1; see Fig. 3). We obtain integrated count 
rates of several to ~ 100 km^^ yr"^ for a depletion factor 
fd — 1/2. Such count rates should be easily detected by 
IceCube, and possibly by AMANDA-II or ANTARES with 
integration times of about a decade (IceCube is planned to 
have replaced AMANDA-II by then). While the character- 
istics of the spectrum presented here are robust, we caution 
that the event rates we obtain are very rough upper lim- 
its, subject to many uncertainties. For example, we have 
assumed that the neutrinos are beamed into the same solid 
angle as the radio beam, which might not be a correct as- 
sumption. The radio beam is thought to be produced within 
~ 107? (see, e.g., Cordes 1978). In our model, the pions are 
produced much closer to the star. They then propagate to 
^ 1000-R before decaying to neutrinos. At this distance from 
the star, the field is not dipolar, and it is difficult to say any- 
thing definite about the distribution of pion trajectories in 
this region. If the neutrinos form a beam, it may be more or 
less collimated than the radio beam. If the neutrino beam is 
more collimated, the neutrino event rates would be higher 
than estimated here. 

The first 807 d of data from AMANDA-II revealed no 
statistically significant sources (Grofi, 2005). Intriguingly, 
there were 10 events (over a background of 5.4) recorded 
from the direction of the Crab pulsar; IceCube will be able to 
confirm or refute this result. While it would be more exciting 
to actually see neutrinos from pulsars, the accumulation of 
null results over the next decade would be interesting as well; 
it would probably mean that photomeson production is in- 
effective or non-existent in the neutron star magnetosphere, 
thus providing a bound on the accelerating potential. 
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